D = 8*ones(1, 1) + 8*eye(3);

L = chol(D)

N = 1000;
x = L * randn(3, N);

mean(x(1)*x(1))/std(x(1, :))
mean(x(1)*x(2))/std(x(1, :))